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Abstract 

^g" ■ The two-body momentum distribution 7/2(^1,^2) of nuclei is studied. First, a compact 

analytical expression is derived for Z = N, ^-closed nuclei, within the context of the inde- 
pendent particle shell model. Application to the light closed-shell nucleus 16 O is included 
and discussed. Next, the effect of dynamical, short-range correlations is investigated in 
the case of the light nucleus 4 He, by including Jastrow-type correlations in the formalism. 
I The effect is significant for large values of p\ and P2 and for angles between the vectors p± 

and f>2 close to 7 =180° and 0°. 

1 Introduction 

The problem of short-range nucleon-nucleon correlations (SRC) in nuclei is still an open 
field in nuclear physics [1]. An indirect evidence for SRC was obtained by the obser- 
vation of a depletion of valence shells in (e,e'p) reactions [2, 3]. It is also well known 
that SRC modify the nucleon momentum distribution [4, 5], i.e. the Fourier transform of 
the one-body density matrix p\{fi,f\i) in the variable r\ — r\>, by introducing significant 
contributions at values of momentum beyond the Fermi momentum. Experimental infor- 
mation from inclusive (e, e') and exclusive (e, e'p) reactions established the existence of a 
high-momentum component in the momentum distribution. 

Beyond the one-body density matrix and momentum distribution, rich information on 
the nuclear ground state and nucleon-nucleon correlations is contained in the two-body 
density matrix (2DM) P2(ri, fy, ry, ry) and its various Fourier transforms in momentum 
space. Experimentally exclusive two-nucleon knockout reactions can provide information 
on the relative behaviour of nucleon pairs embedded in the nuclear medium and give more 
direct insight into SRC [6, 7, 8, 9, 10, 11, 12, 13]. A very important role is played by 
electron-nucleus scattering experiments with the availability of high-energy 100% duty 
cycle electron beams which made possible to carry out double coincidence (e, e'NN) ex- 
periments. Past, present and near future experiments at NIKHEF, MAMI and Jefferson 
Lab use 3 He, 12 C and 16 as targets [8, 9, 10, 11]. 

The analysis of (e,e'2N) experiments is a complicated task [14, 15, 16, 17]. Among 
the most important processes of longitudinal character contributing to electromagnetically 
induced two-nucleon knockout is the scattering off a correlated nucleon pair followed by 
two-nucleon knockout during which each of the two nucleons interacts with the residual 
nucleus (final-state interactions - FSI). The minimal starting point in the analysis is the 
plane- wave impulse approximation (PWIA), in which the FSI are neglected and the two- 
nucleon knock-out cross section can be expressed in terms of the two-hole spectral function 
S(pi ,p2',E). The latter gives the joint probability of removing from the target two nucleons 
with momenta p\ and f>2 leaving the residual nucleus with energy E with respect to 
the target nucleus ground state. Under certain conditions it seems possible to extract 
information on the 2DM and its Fourier transforms and on nucleon-nucleon correlations. 
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Taking the above into account, the development of a simple model of the ground-state 
2DM in order to clarify the two-body characteristics of the nuclear ground state and the 
extent to which it is affected by the SRC appears useful. 

Here we focus on the two-body momentum distribution (2bMD) f]2(pi,P2), which is 
the Fourier transform of the 2DM p2 {r\ , r-i ; fy , ) in the variables f\ — fy and r*2 — ■ 
Another Fourier transform of the half-diagonal two-body density matrix p2hifit ^2, fy) = 
P2(ri,r2;ry,r2) [18], namely the generalized momentum distribution rj(p,Q), has been 
studied previously in the case of infinite nuclear matter [19] and finite nuclei [20]. The 
generalized momentum distribution plays an important role in the study of the FSI. The 
2bMD is straighforwardly connected to the two-nucleon spectral function S(p±,p2; E) via 
an integration with respect to the energy (see Eq.(10) below). The two-nucleon spectral 
function in the case of infinite nuclear matter has been studied using correlated basis func- 
tion perturbation theory including central and tensor correlations [21]. In the case of 16 
it has been evaluated including short-range correlations through defect functions calcu- 
lated via the Bethe-Goldstone equation or correlation functions derived from variational 
calculations and long-range correlations using dressed RPA [15]. 

In this work, first an analytical expression is derived for the 2bMD of Z = N, ^-closed 
nuclei in the independent-particle shell model with harmonic oscillator wave-functions 
and applied to the nucleus 16 0. Our method is an extension of that developed in Refs. 
[20, 22, 23, 24] for the study of the nuclear form factor, the nuclear charge, matter and 
momentum distributions, the one-body density matrix and the generalized momentum 
distribution in closed shell nuclei. This calculation is expected to reproduce the main 
features of the 2bMD, stemming from the finite nuclear size and Fermi statistics, being 
reliable at certain kinematical domains. By introducing the quantity 772 (pi, P2) /vipi) 7 !^) , 
where rj{p) is the nucleon momentum distribution, we pinpoint the effects of the finite size 
and of the statistical correlations. 

Next, the effect of the SRC on the 2bMD of 4 He is studied by including Jastrow-type 
correlations in our calculation, via the lowest term of a cluster expansion. This so-called 
low order approximation (LOA) [25] has been exploited for the one-body density matrix 
and the two-body density matrix by Bohigas and Stringari [26] and Dal Ri, Stringari 
and Bohigas [27] and has been widely used to study single-particle nuclear properties, 
such as the (point-nucleon or charge) density [28, 29], the elastic or charge form factor 
[26, 27, 29], the one-body density matrix in coordinate [28, 30] or momentum space [28] and 
the nucleon momentum distribution [26, 27, 28, 30], as well as two-body nuclear quantities, 
namely the two-body density matrix in coordinate space [31], the static structure function 
S(Q) [27], the relative pair density distribution ^2(^12) [27, 31], the center-of-mass pair 
density distribution P2{Rcm) [31], the two-body relative momentum distribution 772 (& e i) 
[31, 32], the two-body center-of-mass momentum distribution ^{Pcm) [31, 32] and the 
combined two-body center-of-mass and relative momentum distribution i]2(Pcm, <7rci) [32]. 
A two-gaussian central correlation function, f(r) = 1 — c\ exp(— r 2 /(3 2 ) + C2 exp(— r 2 //3|), 
is used in this work. We find significant deviations from the independent-particle picture 
at high momenta, especially when the vectors p\ and f>2 are antiparallel or parallel, as it 
is clearly shown by the values of the quantity V2(pi,P2)/v(Pi)v{P2)- There is, however, 
within the LOA, significant sensitivity of the results in some kinematical regions, mainly 
on the correlation function used to describe the SRC. For certain correlation functions one 
may obtain negative values of 7/2 in few regions of momenta. 

The definition and properties of the two-body momentum distribution are presented 
in detail in Section 2. In Section 3 an exact expression is derived for the 2bMD of Z = N, 
^-closed nuclei in the independent-particle shell model with harmonic oscillator wavefunc- 
tions and applied to the case of 16 0. In Section 4 the effect of SRC on the 2bMD of 4 He is 
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explored by including Jastrow-type correlations in our calculation via the LOA, and the 
results as well as the reliability of our approximation are discussed. Finally, in Section 5 
we give a summary of the results and hints for possible future developments. 



2 The Two-Body Momentum Distribution 

The two-body momentum distribution 772(^1,^2) gives the combined probability density 
of finding in the nucleus a nucleon with momentum p\ and another one with momentum 
P2- In a system of A (A > 2) identical particles, described by a unit-normalized state |^), 
it can be defined as the expectation value 

m(Pi,P2) = (^T.^sAu^^P^)- (!) 



s.s' 



The 2bMD is the diagonal element of the 2DM in momentum space V2(pi,P2',Pi' ,i>2'), 

m(pi,P2) = m{pi,P2-,pi,P2) (2) 

and may also be defined as the Fourier transform of the 2DM p2 {f\ , r*2 ; fy , ) in the 
variables r\ — ry and r*2 — '■ 



m 



(pi,P2) = J ^(^i.^jri/.faOe-^-^-^Je-^-^-^^dVidVadVdVa/. (3) 



(By taking h = 1, the 2bMD has the dimension (length) 6 .) The normalization adopted is 
such that 

/ %(pi,p 2 )d 3 pid 3 p 2 = A(A - 1), (4) 
and the following sequential relation holds: 

J m(Pi,P2)d 3 p 2 = (A - l)r/(pi), (5) 

where rj(p) is the one-body momentum distribution, normalized to the number of particles 

J 7](p)d 3 p = A. (6) 

Equation (5) connects the 2bMD to the generalized momentum distribution [19, 20] 
(V(P,Q = 0) = (A-I) v (p)). 

The 2bMD 7/2(^1,^2) is related to the two-body center-of-mass momentum distribution 
V2(Pcm) (Pcm = Pi + P2 is the two- nucleon center-of-mass momentum), and to the two- 
body relative momentum distribution 772 (^rei) (ft-ei = \{pi —P2) is the relative two-nucleon 
momentum) [31, 33] 



m 



{Pcm) = j m{PcM/2 + q ieh P CM /2 - g rc i)d 3 g re i, (7) 

(ftcl) = j m(PcM/2 + q rch Pcm/2 - fei)d 3 P C M- (8) 

It is also related to the combined two-body center-of-mass and relative momentum distri- 
bution rj 2 (Pcm, qvei) [32] 

m( p CM, <?rei) = J m(PcM/2 + q Teh P C m/2 - g rc i)dQp CM dfi 9rcl . (9) 
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Its relationship with the two-nucleon spectral function S(pi,p2", E) is 



m(pi,P2) = J dEStfufcE). (10) 
In a totally uncorrelated system the 2bMD is factorized as 

m(Pi,P2) = (i - -^MpiMrO- (11) 

The constant 1 — 1/A appears due to the normalization condition (4). For the ground 
state of a non-interacting Fermi system (a system of A identical, non-interacting fermions) 
where only statistical correlations are present, starting from the 2DM [4], the 2bMD can 
be shown to have the form 

m(pi,P2) = v(pi)v(p2) - -\m(pi,P2)\ 2 , (12) 

where r}i(pi,P2) is the one-body density matrix in momentum space and v is the level 
degeneracy. Note that in the special case of a non-interacting Fermi system with A = v 
Eq. (12) becomes identical with Eq. (11). 

As an example we consider the ground state of an infinitely extended ideal Fermi gas 
(volume Q and number of particles A going to infinity while the particle density stays 
constant) with Fermi wave number kp and degeneracy v. Equation (12) yields for the 
2bMD per particle-pair 7% (pi,P2) = V2 (Pi->P2)/[A(A — 1)] 

v!(Pi,P2) = 2 0(k F - Pl )e(k F - P2 )[l - 5 PlP Jv\. (13) 

The momentum distribution per particle, if \p) = r] F (p)/A = r] F (p)/A, is given by 

f(p) = (^*4) ^(kp-p). (14) 

Let us define the quantity £ as a measure of finite-size effects, of statistical and short-range 
correlations of dynamical origin 

t(pi,P2) = m(pi,P2)/v(pi)v(p2)- (15) 

For an infinitely extended ideal Fermi gas £ is defined for pi,p2 < kF and equals 

t F ( -> - n f 1, Pi ¥"P2 

? (Pl,P2) = < , , I - - • (16) 

y 1 - l/u, pi=p 2 

From Eqs. (11) and (16) we realize that, in the case of the infinitely extended ideal Fermi 
system, if p\ 7^ P2, £ equals 1. This holds even in the presence of dynamical SRC, if 
long-range order does not exist [34]. In the case of the infinitely extended ideal Fermi gas, 
£ = 1 — l/u \i p\ = §2- As for the finite, non-interacting Fermi system, we realize from 
Eq. (12) that, if p\ = P2, £ = 1 — 1/za Deviations of £ from this value show the effect of 
dynamical correlations. For p\ 7^ P2 deviations of £ from 1 — 1/^4 is a measure of statistical 
and (or) dynamical correlations in a system of finite size. 
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3 Method of Calculation for ^-Closed Nuclei - Ap- 
plication to 16 

3.1 Method of Calculation 

Let us now consider a system of A identical non-interacting fermions in its ground state. 
The fermions occupy the lowest single-particle energy eigenstates \rij) (j = 1, 2, Ajv) 
and the 2bMD is given by Eq. (12). In the case of the nucleus, A is the mass number 
and v is the degeneracy due to the nucleon spin and isospin.We consider Z = N, ^-closed 
nuclei. In order to obtain closed analytical expressions for rj2(pi,P2), we have assumed 
that the nucleons move in an isotropic harmonic oscillator potential and that the spin-orbit 
coupling is negligible, the center-of-mass and finite nucleon size corrections are small and 
the Coulomb interaction (relevant for protons) is small. Then Eq. (12) gives: 

ax 

m(Pi,P2) = ^e-^V^ 2 J2 (Pi 6 )" 1 E (P2&r<Wcos 7 ), (17) 

Ml=0 A*2=0 

where the coefficients G^^ are given by 

Ml M2 

Gmna ( cos 7) = 4/^/2/^/2 ~ E E **iAa ( COS T)^( M i-Ai)( M2 -A 2 ) ( cos 7) (18) 

Ai=0A 2 =0 

and the coefficients / M / 2 and /Ca 1 a 2 ( cos 7) are given by 

U/2 = E /nf> (19) 
n£,occ. 

/C AlA2 (cos 7 )= £ ^ lA2 ^(cos 7 ). (20) 

n£,occ. 

In the above equations 7 is the angle between the vectors p\ and P2, b is the harmonic os- 
cillator parameter, Pt(x) is a Legendre polynomial and iV max = (2n + £) max is the number 
of energy quanta of the highest occupied n£-level. The coefficients /^/ 2 , K^ 2 are ratio- 
nal numbers that enter the corresponding expressions of rj(p) and r]i(pi,py) respectively. 
Analytical expressions for f^ 2 and K^ 2 can be found in Ref. [20]. They are different 
from zero only if the indices \x and Ai + A2 are even. 

The corresponding expression for the spherically symmetric nucleon momentum dis- 
tribution rj(p) is [22] 

V(P) = v(p) = -^ e ~ p2&2 E (pbf'Vx- (21) 

71 A=0 

It is verified that the 2bMD 172(^1, P2) given by Eq. (17) satisfies the property of Eq. (5). 

3.2 Application to ie O 

Using the expressions (17) and (21), we obtain results for the quantities rj2(pi,P2) and 
£(pi,P2), in the case of the nucleus 16 O. The harmonic oscillator parameter b = 1.7825fm 
has been determined so as to reproduce the experimental value of the charge root mean 
square radius of 16 0, (r) 1 / 2 = 2.737fm [35]. In Fig. 1 the 2bMD of 16 in the harmonic 
oscillator model (continuous lines) is plotted for p\ parallel to P2 and for p\ = 0, 1 and 
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1.5 fm _1 , as a function of the variable P2 (P2 = P2P1). In the same graphs the 2bMD of 
an infinite ideal Fermi gas is shown (dashed lines), calculated via Eq. (13) with v = 4 and 
with Fermi wave number kp equal to the one corresponding to 16 0, namely kp = l.lfm -1 . 
This latter value has been calculated using local density approximation and coincides with 
the one used in Ref. [36] . In order to make the comparison meaningful, the 2bMD of the 
ideal Fermi gas was normalized to the same number of particle pairs as those found in 
16 0, namely to A(A - 1) = 240. 

It is evident from Eq. (13) that the 2bMD of the ideal Fermi gas exhibits disconti- 
nuities at pi 5 2 = kp and at p\ = Pi- This is observed in the first two panels of Fig. 1 
{pi = 0, 1 fm" 1 ) at P2 = pi and at P2 = ±kp. In the case of the finite system 16 O in the 
harmonic oscillator model a rather similar behaviour is observed, but the 2bMD changes 
smoothly with P2. For p\ = 1.5fm _1 (which is higher than kp = l.lfm" 1 ) the 2bMD of 
the ideal Fermi gas has vanished and that of 16 O has dropped considerably. 

In Fig. 2 the quantity £(pi,P2) is plotted in the case of 16 O as a function of cos 7 and 
for pi = P2 = 0, 1 and 4 fm _1 . The values of £ lie between 0.75(= 1 — and 1. £ is 
equal to 0.75 when p\ = f>2 (as shown for example in the first panel, and in the other two 
panels when cos 7 = 1). We realize that for high values of p\ and P2, £ tends to 1 and 0.75 
for cos 7 = and cos 7 = —1 respectively. The deviation of the values of £ from 1 or 0.75 
is a finite-size effect in the presence of statistical correlations. In the case when p\ 7^ P2 
(not shown) the quantity £ as a function of cos 7 remains higher than 0.75 and below or 
equal to 1, following a similar behaviour as in Fig. 2, namely £(—1) > £(1) and £(#o) = 1 
for some value cos 7 = xq < 0. 



Up to now we have ignored the effect of dynamical, short-range correlations on the 2bMD. 
It is expected that such a calculation is reliable only in restricted kinematical domains, 
namely for p±, P2 < kp. In this Section we make an attempt to investigate the effect of 
the SRC by including Jastrow-type correlations in our study. 



The approach adopted here is based on the Jastrow formalism with a state-independent 
central correlation function, which introduces the short range correlations. Employing 
the low-order approximation (LOA) of Ref. [27] for the 2DM and performing the spin- 
isospin summation (y = 4) in Eq. (14) of Ref. [27], we obtain for the correlated 2DM 
P2° vv (fi, r*2; fy, r 2 ') the following expression: 



where r^j = \fi — fj\, g(r,r') = f(r)f(r') — 1 and f(r) is the correlation function, which 
has to obey the conditions /(0) < 1 and for r — > 00, f(r) — > 1. 



4 The Effect of Short-Range Correlations 



4.1 Method of Calculation 




(22) 
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In Eq. (22) the one-body density matrix pi(r,r'), the 2DM P2(ri, rg; fy, ry) and the 
half-diagonal 2DM p2hif\i'?ii ^i') are calculated in the harmonic-oscillator model. The cor- 
related 2bMD, r]2° rr (pi,P2), is then calculated by Fourier transforming P2° rr (^i>^2;rL',^2') 
according to Eq. (3). The correlated momentum distribution r] corr (p) is calculated likewise 
in the LOA by performing the spin-isospin summation in expression (13) of Ref. [27] for 
the one-body density matrix and by Fourier-transforming with respect to f\ — fy. 

The LOA preserves the normalization of the density matrices. As a consequence, the 
one- and two-body momentum distributions calculated in the LOA obey the normaliza- 
tion conditions (6) and (4), as well as the sequential relation (5). However, in the LOA 
some probability distributions (i.e. positive-definite quantities) may obtain negative values 
within regions of their domain. This problem is probably related to the A— representability 
of the LOA two-body density matrices. For example, it was encountered in calculations of 
the relative distribution function p 2 (ri 2 ) [27] and the combined two-body center-of-mass 
and relative momentum distribution ^{Pqm-, Orel) [32], which implies that it may show up 
in the 2bMD too, since the latter is related to 772 (Pcm, q Te i) via Eq. (9). In the case of the 
^2(^12) the problem was subsequently solved by an appropriate modification of the corre- 
lation function [31]. The above remarks are highly relevant to the choice of the correlation 
function used in this work, as we will discuss in due course. 



4.2 Application to 4 He 

The 2bMD and the momentum distribution of 4 He in the harmonic oscillator model, 
which we will hereafter denote by r/ 2 (pi , P2 ; b) and rf (p; b) respectively, are given by (see 
Eqs. (17), (21)) 

126 6 

7T 3 



V^Pi,P2-M = ^~ b2[pl+pl \ (23) 



^ 6 ) = ^ e ~ p2b2 - ( 24 ) 

yielding £(pi,p~2) = 0.75. The case falls into the category of non-interacting Fermi systems 
with A = v, hence Eq. (11) holds (see remark under Eq. (12)). The harmonic oscillator 
parameter b = 1.382fm reproduces the experimental value of the charge root mean square 
radius of 4 He, (r 2 ) 1 / 2 = 1.67fm [35] in this model (corrections due to center-of-mass motion 
and finite nucleon size are taken into account). We include SRC in the LOA using a two- 
gaussian (2G) correlation function, f(r) = 1 — c\ exp(— r 2 / [3 2 ) + C2 exp(— r 2 //3|). The 
expression of the 2bMD is 

V2° rr (Pi,P2) = V2(Pi,P2]b) +6rj 2 {pi,P2;b,c 1 ,y 1 ) + 5r] 2 (pi,p 2 ; b, -c 2 ,y 2 ) 

_ 24c lC2 -| - e-^^+^tf^+^+i&i***] 

1 2 vr 3l [(l + 4 yi )(l+4y 2 )] 3 / 2 

2 l+2 yi + 2y 2 2 h 2 2 ,2 

+ ?m 1 q W2"[ e 1+3 ^ +3 ^+^l^ Pl e" P 2 fe +(P1^P 2 )] 

(1 + 3yi + 3y 2 + SymY 12 

e -fe 2 ( P f+p|)| ( 25 ) 



(1 + 2j/i + 2y 2 ) 3 / 2 
where yi = b 2 /(3 2 and 

r /- - 1. \ 12 ^ 6 rr 10c ^C 2 _ b 2, 2, 2^ 

S m ( Pl , P 2;b,c,y) = —{\jYTWT 2 ~ ^TW^ ]e 



2C »,2r l+3y / 2 I 2\ I 



2 II 



(1 + 4y) 3 / 2 
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+ 



(1 + %)3/ 2 



4c 



.[ e -TTip?^ 2 e -pl& 2 + fa ^ P2 )]}. (26) 



(l + 3y) 3 / 2 

The corresponding expression for the momentum distribution r/ corr (p) of 4 He is 
?7 corr (p) = ?? (p;6) +*r/(p;6,ci,yi) + Srj(p; b, -c 2 , y 2 ) 

„2 h 2 l+2K 1 +2y 2 



-6cic 2 



7T 



3/2 



where 



5r)(p;b,c,y) 



1 + 3yi + 3y 2 + 8 yi y 2 ) 3 / 2 (1 + 2 yi + 2y 2 ) 3 / 2 
126 3 



7T 



3/2 



ere 1+2 f 2ce 1 + 3 f 

[(l + 2y)(l + 4y)3/2 " (1 + 3y) 3 / 2 



(27) 



(28) 



If we put c 2 = in the above equations we end up with the corresponding expressions for 
the single gaussian (1G) correlation function, f(r) = 1 — cexp (— r 2 / 1 ) which has been 
used in evaluating several quantities of 4 He [26, 27, 28, 29, 31, 32]. We have introduced 
the 2G correlation function, that can overshoot unity at intermediate distances. This was 
found necessary in order to eliminate the negative values of ?? 2 (pi,p 2 ) that result at some 
regions of intermediate values of momenta whenever a 1G correlation function is used. 

The parameterization of the new, 2G correlation function by means of a rigorous pro- 
cedure such as minimization and fitting to some set of experimental data (e.g. charge form 
factor data) lies beyond the scope of this work, therefore it was obtained somehow heuris- 
tically. First, it has been chosen to maintain the small-r behaviour of the 1G correlation 
function of Ref. [31], by setting c\ — c 2 = 0.8 and [3\ = 0.8fm. Subsequently, the extra 
parameters have been determined requiring that the healing distance is equal to or some- 
what lower than lfm and the relative pair density distribution p 2 (ri 2 ) is approximately 
zero at ri 2 = 0. The value of the harmonic oscillator parameter b was again chosen so as 
to reproduce the experimental charge root mean square radius of 4 He. We ended up with 
c 2 = 1.7 (and therefore c\ = 0.9), /3 2 = l.lfm and b = 1.4fm. In the left panel of Fig. 3 the 
resulting 2G correlation function, that overshoots unity, is plotted along with the 1G one of 
Ref. [31]. The corresponding values of a wound parameter k = J(l — /(ri 2 )) 2 /C>2(Ti2)d 3 ri 2 
(where /9 2 (ri 2 ) is the relative pair density distribution function calculated in the harmonic 
oscillator model and normalized to unity) for the 2G and 1G correlations functions are 
equal to 0.008 and 0.018 respectively. In the center panel of the same figure the momen- 
tum distribution per particle rj{p)/A is plotted for both correlation functions, 2G and 1G, 
along with the result of a variational calculation [36, 37]. In the right panel of Fig. 3 the 
proton-proton relative pair density distribution (normalized to unity) p pp {r 12) is shown 
for both correlation functions 2G, 1G. The bars cover the results of several variational 
calculations using a variety of two- and three- nucleon interaction models [38]. The values 
of rj(p) and p pp (ru) calculated with the two-gaussian correlation function compare quite 
well with the ones of other methods. 

Some numerical results for ?? 2 (pi,p 2 ) are presented in Figs. 4 and 5. In Fig. 4, the 
2bMD of 4 He is plotted in logarithmic scale, for p 2 parallel to p±, p 2 =p 2 pi, as a function 
of P2, for p\ = 0, 1, 2, 3 fm _1 . The correlated 2bMD, as given by Eq. (25), is plotted 
with continuous lines and the harmonic-oscillator result, Eq. (23), with dashed lines. It 
seems that, within the present approximation, the correlated 2bMD shows deviations from 
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the harmonic-oscillator picture mainly at high values of p\ and/or P2, developing high- 
momentum tails. In most cases, the deviations are larger when pi, pi are antiparallel 
rather than parallel. 

In Fig. 5 the quantity £(pi,P2) f° r 4 He is plotted in logarithmic scale, as a function 
of cos 7 for pi = p2 = 2 and 4 fm _1 . The continuous line includes SRC and the dashed 
line corresponds to the harmonic oscillator model and it equals 0.75. Deviations from this 
value shows the effect of SRC. The correlated £ in the present approximation grows larger 
than unity at backward angles, and more so for large values of pi, p2, as indicated by 
the pi = p2 = 4ml" 1 graph. In this latter graph, deviations from 0.75 are large even at 
forward angles. We should note that for low values of p\ and p2 (not shown), the correlated 
£ takes values between about 0.7 and 1. 

We have also compared our results for the 2bMD with those obtained by using the 
single-gaussian correlation function of Ref. [31]. Although the overall behaviour of the 
2bMD is not dramatically affected by the choice of f(r), the two sets of results differ 
largely at certain kinematical regions, in particular at intermediate values of momenta 
where the single-gaussian correlation function yields negative results. (If the 2DM is ap- 
proximated by only the first terms of the expression (22), i.e. p™ 11 = /(^i2)/( r i'2')P2) 
these regions of negative values show up as exact zeros at certain points [39].) As antic- 
ipated in Sec. 4.1, a similar, unphysical situation has been encountered in the study of 
other two-body quantities, when the LOA is used [27, 32]. It has been argued [32] that 
such regions indicate where the SRC start to show up dramatically, and are most sensitive 
to the approximation employed to account for the SRC. The results in these regions are 
also sensitive to the particular correlation function used. Let us note here that even in the 
case of the 2G correlation function a somewhat lower value of the oscillator parameter (eg. 
b = 1.35fm) retrieves the negative values of the 2bMD, indicating that the uncorrelated 
single-particle wavefunctions also play some role. 

5 Summary and Conclusions 

In summary, a compact analytical expression has been derived for the two-body momentum 
distribution, r]2(pi,P2), of Z = N, ^-closed nuclei within the context of the independent- 
particle shell model, using harmonic-oscillator wavefunctions. An application to the 16 
nucleus revealed interesting features arising from the Fermi statistics and the finite nuclear 
size. Next, the effect of short-range correlations was investigated by including Jastrow- 
type correlations in the calculation of T)2(j>i,P2) of the 4 He nucleus, using the low-order 
approximation of Ref. [27] and a two-gaussian correlation function. Significant deviations 
from the independent-particle picture were found for large values of p\ and P2 and for angles 
between the vectors p\ and P2 close to 7 =180° and 0°. The quality of our approximation 
was also explored. It was found that, at certain regions of momenta, the results are 
sensitive to the correlation function used. 

Using the same approximation, one could evaluate the 2bMD and investigate the ef- 
fects of short-range correlations in other interesting nuclei, for example in 16 and 40 Ca. 
Furthermore, a general expression within the harmonic oscillator model applicable to ev- 
ery I— closed nucleus could be derived for the two-body momentum distribution 772(^1,^2)) 
the two-body center-of-mass momentum distribution 772 (Pcm) and the two-body relative 
momentum distribution r/2(g rc i), as well as for the two-hole spectral function S(pi,p2',E). 
The above calculations are relevant to the study of the two-nucleon knock-out reactions. 
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Figure 1: Two-body momentum distribution for p\ parallel to f>2, as a function of p 2 
(p 2 =P2pi) for pi = 0, 1, 1.5 fm _1 , in the case of 16 in the harmonic oscillator model 
(continuous line) and of an infinite ideal Fermi gas with k-p =l.lfm _1 , normalized 
to the same number of particle pairs as in 16 O (dashed line). 




Figure 2: The quantity £(Pi>P2) (see Eq. (15)) as a function of cos 7 (7 is the angle 
between the vectors p\ and p 2 ) for Pi = P2 = 0, 1, 4 fm _1 , in the case of 16 in the 
harmonic oscillator model (continuous line). 
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Figure 3: Left: The single-gaussian correlation function with c = 0.76, j3 = 0.813fm, 
from Ref. [31] (1G CF) and the two-gaussian correlation function employed in this 
work with the parameterization described in the text (2G CF). Center: The mo- 
mentum distribution per particle of 4 He calculated in LOA with the 1G and 2G cor- 
relation functions is compared with the variational calculation of Ref. [36] ("Var."). 
Right: The proton-proton relative pair density distribution (normalized to 1) of 4 He 
calculated in LOA with the 1G and 2G correlation functions is compared with the 
domain of results obtained via a set of variational calculations [38] ("Var."). 
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Figure 4: Two-body momentum distribution (in logarithmic scale) for p\ parallel 
to p 2 , as a function of p 2 (j>2 —ViPi) f° r Pi — 0, 1> 2, 3 fm _1 , in the case of 4 He 
including SRC, Eq. (25) (continuous line), and in the harmonic oscillator model, 
Eq. (23) (dashed line). 




Figure 5: The quantity £(pi,P2) (see Eq. (15)) in logarithmic scale, as a function of 
cos 7 (7 is the angle between the vectors p\ and p 2 ) for pi = p 2 = 2, 4 fm _1 , in the 
case of 4 He including SRC (continuous line), and in the harmonic oscillator model 
(dashed line). 
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